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1. 

Underwater  Explosion  Bubbles  II: 
The  t^ffect  of  Gravity  and  the  Change  of  Shape 
Ignace  I.  Kolodner  and  Joseph  B.  Keller 

I    Introduction 

The  classical  theory  of  an  underwater  explosion 
bubble  is  based  on  the  assumption  that  the  bubble  re- 
niains  spherical  at  all  tln=?s .   However  -actual  bubbles 
do  not  remain  spherical  although  they  rray  be  spherical 
Initially.   Instead  they  become  flattened  or  even 
kidney-shaped  and  often  break  up.   This  change  of 
shape  is  mainly  due  to  gravity.   The  present  report 
presents  a  theory  which  describes  the  change  of  shape 
as  well  as  the  rise  of  the  bubble,  by  taking  account 
of  grsvlty. 

Several  other  authors,  such  as  H  rring  and  Ward, 
have  also  riven  theories  of  the  chpnge  of  bubble  shape 
^u"-  t'""  frravlty.   The  ■prese'^t  worT-r  is  Int'^'nded  to  be 
mere  system.atlc  and  more  complete  than  any  of  the 
former  theories. 

In  the  present  theory  it  is  assumed  that  the 
writer  is  incompressible  and  unbounded,  and  the  pressure 
is  assumed  to  be  uniform  throughout  the  bubble  at  all 
times.   Initially  the  bubble  is  assuried  to  be  spherical. 
Then  if  we  neglect  gravity  wo  find  that  the  bubble 
remains  spherical  but  performs  periodic  radial  oscil- 
lations, exactly  as  in  the  classical  theory.   Taking 
account  of  first  order  terms  due  to  gravity,  we  find  that 
the  bubble  remains  spherical  but  also  rises  exactly  as 
predicted  b7/  Herring's  formula.   Higher  order  terms 
lead  to  the  change  of  shape  and  sometimes  to  the  break- 
up of  the  bubble,  and  also  lenH  to  a  mcdlflc^tlon  of 
the  Herrinn  rise  formula.   ^11  of  these  results  are 
discussed  in  section  VIII  of  this  report. 

The  assumption  that  the  water  is  incompressible 
can  be  removed,  as  it  was  in  our  previous  report  on 


iori-l 


9.. 
the  spherical  bubble.   However  it  does  net  seem  worth- 
while to  add  the  cc-^iipiic-.ticn  lOf  corrprossibil  ity  to  the 
other  difficulties  of  the  present  problem,   f^urthermcre 
in  section  VIII  v/e  have  indicated  how  the  main  effects  of 
compressibility  can  be  obtained  by  combining  the  present 
results  with  those  of  our  previous  report. 

The  assumption  th^t  the  wf^ter  is  unbounded  cnn 
also  be  removed,  and  the  theory  without  this  assumption 
is  the  subject  of  our  next  report.   In  that  report  it 
will  be  shown  that  the  results  of  the  present  report 
provide  a  first  apprcximaticr  to  the  solution  when 
the  bubble  is  nf-^t  too  near  the  bcund^ries  .   Of  course 
further  approximations  will  also  '-_:  rriven. 

In  section  II  the  prob]em  is  f  .-^rmul^ted ,  in  section 
III  mxvinR  coordinates  are  introduced  and  in  section  IV 
dimensicnless  variables  are  defined.   The  method  of 
S'^lution  is  explained  in  section  V  and  the  solution  up 
to  and  includine  the  thir'-i  order  is  considered  in 
section  VI.   A  special  c^^se  is  treated  exnlicitly  in 
secticn  VII.   Finally  in  section  VIII  all  the  results  are 
discussed  and  some  general  conclusions  are  drawn. 

II   Fcrm.ulaticn 

We  assume  that  an  incompressible  invistid  fluid 
of  infinite  extent  crntrins  a  fras  bubble  within  it. 
The  velocity  u(x,y,'^,t)  of  the  fluid  is  assumed  to  be 
derivFible  from  h  notenti^l  function  (|)(x,y,z,t)  which 
satisfies  Laplace's  enuati'~n 

(1)  u^=V(i> 

(2)  ^^4=   0 

The    pressure    p(x,y,-,t)in   the   w^ter    is    then   (riven  b:/   the 
Bernoulli    enuation 


(3)       P   =  m(t)    -f  Ut    -^  f(V4)^_ 


f^- 


pz 


r:   •    .     I 
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3. 

In  (3)  P  is  the  density  of  water,  g  is  the  acceleration 
of  gravity  and  m(t)  is  an  arbitrary  function  of  time. 
We  have  also  assumed  that  the  positive  z  axis  points 
vertically  upward. 

The  velocity  V^    is  assumed  to  vanish  at  infinite 
distance  from  the  bubble 


[k)       V<P  ^0   as  (x  +  y  +  z  )  ' 


CO 


From  (1)  it  is  clear  that  ^    is  defined  up  to  an 
additive  function  of  t.   This  indefiniteness  can  be 
removed  by  specifying  ^   or  (j).  at  one  point  for  all  t. 
For  this  purpose  we  assume 

P    P    9  1  /"^ 

(5)  (!>-*►  0    as(x  +  y'+z)'   -->0D  . 

The  condition  (5)  implies  that  (l^^-^^O  at  infinity,  and 
also  im.plies  (4),  which  can  therefore  be  omitted. 

The  arbitrary  function  m(t)  in  (3)  can  be  deter- 
mined by  specifying  p  at  ono  point  for  ^11   t.   For 
this  purpose  we  assume  that 

2    2 

(6)  p(x,y,o,t)->  P^   as  (x  +  y  ) ->aD 

Thus   m(t)    can  be   determ.ined    from    (3,5,6)    vrhich   yield 

(7)  m(t)  =  P^. 

We  call  P   the  pressure  at  infinity  at  z=0.   [It  is 

o 
just  the  hydrostatic  pressure  at  the  level  z=0, 

namely  P  =  p  +  f gz  ,  if  p   is  the  atmospheric  pressure 
"'oo-'o      o 

and  z  >o  is  the  depth  of  z=o  below  the  water  surface.] 
o 

The  bubble  surface  is  assumed  to  be  given  by  the 
equation 

(8)  F(x,y,z,t)  =  0 

Since  the  normnl  velocity  of  the  bubble  surface  m_ust 
be  the  same  as  that  of  the  nd.i'acont  water,  F  must 
satisfy  the  kinempttic  condition 


1  ("1  .■       V  ,         I 
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(9)  V^F  •  Vcj)  +  F   =  0    on   F  =  0. 

The  pressure  it  within  the  bubble  is  assumed  to 
be  a  known  function  of  V,  the  bubble  volume.   'Ve  assume 
that  this  function  is  the  adianatic  one,   Tt  =  KV'  , 
where  K  and  5  are  constants,  5  bein^  the  adiabatic 
exponent  for  the  sas  within  the  bubble.   Then  because 
pressure  must  be  continuous  across  the  bubble  boundary, 
we  have  the  dynamic  condition  p  =  KV    or,  usine  (3) 

(10)  P^  -  ^h+     ^    {^^)j       -   ^gz  =  KV'^   on  F  =  0. 

The  msthematical  problem,  which  we  consider  is 
that  of  finding  ^    and  F,  satisfying  (2,5,9,10),  given 
F(x,y,z,o)  and  F  (x,y,z,o).   It  is  not  necessary  to 
specify  ^   initially  since  it  is  determined  by  (2,5,9)* 
In  particulfir,  we  will  assume  th?t  the  bubble  is  ini- 
tially 8  sphere  of  raduis  A   and  that  its  initial 
radlpl  velocity  is  a  constant  ^  .   Once  i  is  deter- 
m-lned ,  u   and  p  are  given  by  (1,3). 


Ill   Coordinates 

Suppos°  that  the  ori.crln  of  the  x,y,z  cccrdinate 
system  is  fixed  at  the  center  cf  the  bubble  at  t  =  0. 
We  nov/  introduce  a  moving  coordinate  system  f,V^,2;  by 
the  eouations 

(11)    f=x,V-^=y,  ?'  =  z-B(t). 

In  (11)  B(t)  is  a  function  to  be  determ.ined  subseauent- 
ly  in  such  a  way  that  the  orisin  is  always  at  the 
center  of  gravity  of  the  bubble.   Kevins  c" "rdinftes 
are  introduced  in  order  that  the  criccin  remain  inside 
the  bubble  :  s" I'ng  as  possible. 

Because  of  the  cylindrical  symm.etry  of  the 
problem  it  is  advantageous  to  introduce  the  spherical 
coordinates  r,0,co.   Obviously  the  solution  will  be 
independent  of  co ,  and  we  will  assume  this  from  the 
outset.   We  define  r,9,co  by 


(        1 


(12)  r   =  r+    ^     +   r   =  x^+   y  •  +  (z-B)"    ,    r>0 


COS0    =    ■£    =     I  r^-. ^7~~7X^  .        Oi©<7t 


£  =  .A-p — vrr-^ 

r    V  X  +  y  +(z- 


y  '+(z-B)' 


» 


f     X 

tanco  =  -rr  =  —  ,   0<a3<2iT 

\       J 

Now  Introducing  R(©,t)  we  write  the  eauation  of  the 
bubble  surface  in  the  form 

(13)    F(x,y,z,,t)  =  r-R(e,t)  =  0 
The  bubble  volume  is  then  given  by 

Uk)  V(t)  =  ^      r'^   R^(©,t)sin9d0 

^o 
The  condition  that  the  moving  origin  be  at  the  center 
of  prravity  of  the  bubble  may  be  written  as 
r 


(15)  {       R'^cos9d9  =  0 
■^o 

To  express  (P,9,10)  in  terms  of  r  and  9,  we 

observe  thf^t  if  l(£,^,^,t)  is  the  expression  in  terms 

of  E,^,i::   for  seme  function  I(x,y,z,t),  i.e.  if 

I(x,y,z,t)  =  T(e,>^,^:,t),  then  wo  have  the  foil-wing 

differentiation  formulf.e: 

(16)  1^=1^,  Iy=I.,  I^=T^,  I,=I,-RI^-I^-B(cr.s9T^.s|n9Tg) 
If  we  apnly  (16)  to  (9,10)  we  obtain 

(17)  k'^-4.   *9  =  R^+ft(cos9+^^^l^Rg)    St  r  =  R(9,t) 

(18)  K(|l  5^  R^\in9d9)-^  ^  P^^-J^|i(rcos9+R)Ht-^(cos9^^. 

sinQ       -   p  -    ~ 
-r    (j)  )+  ^(^^-k:_(J)  )l   at  r=R(9,t). 
w    ■-   ^  r    ^  ;' 


e'r-*.^        ^r,      -,   ., 


■-.  ;  :~ 


),;,-,  ^U 


6. 
The  problem  is  now  that  of  fin'iine-  (f),R{0,t) 
and  B(t)  satisfying  (P,5 ,15 , 17 , 18)  and  the  initial 
conditions  (19) 

(19)  R(9,^)='^o,  Rt(e,o)='-^,  B(g)=B(o)=0 

It  will  prove  convenient  to  employ  the  energy 
equation,  which  Is  a  direct  consequence  of  the  pre- 
ceding equations,  namely 

(20)^^(^  C   R^sinGdG)^"^  +  ^(  P  -rgB)  cVsinOdG  - 
-1      >c  --c 

-ttK  H'^(Rt+B[cos94-^|^  RQ](|)(R,9,t)sin9d©  = 
'  o 

The  terms  on  the  left  side  are  the  internal  enerp-y  of 
the  bubble,  the  potential  energy  and  the  kinetic  energy; 
of  the  water,  respectively:  the  terms  on  the  risht  are 
the  correspondincr  auantlties  evaluated  at  t  =  0,  and 
E  is  the  total  enersy  which  is  defined  by  this  eouaticn. 


IV   Dimensi  nless  Variables 

It  is  crnveni3nt  to  define  dimensicnless  vari- 
ables r,'t,f,X,b,r',a  ,0^'-   by  introducing  a  unit  of 
length  L  and  a  unit  of  time  T.   These  quantities  are 
define'-^  as  follows: 

{21)r=L7      B(t)  -  Lb(t)     k  =  ^V   " 


5-1 
■Kr.e.t)  =  L''T"^i;(?,e,t)   A  =  La  "     k=  -^-(^zl)* 

^  GO  _-. 


6-1   6 


R(0,t)  =  LX(©,t)    l^=   LT"^a 


O  ■' 


f       : 


7. 
For  L  and  T  we  choose  the  definitions 

Now  equations(2,5,l!^>17,l8,19)  become  the  following 
set  for  "^ ,\   and  b: 

(23)  V  4  =0j   (j>->Oasr— >oo 

(2i|.)g(e,t)  =  4^  -  -|  4^  -  X^  -  b(cos  9  -:-  -jr-   >^q)  "  0 

r 

at  r  =  X(9,t) 
(25)h(9,t)  f  4t  ^-  I  ([1  +(-x)^H  ^  4q  +  b  sin9]2  + 

+2X^  ^  [  i  4g  -'-  b  sin9]  +  X^  -  b^y+  <^\?^+   b)  + 


^(i  (   X^sinGd©)    -1=0   at  r  =  A(9,t). 
kl2  ..L         J 


+   -(- 


■^  o 


y 


(26)  X(9,o)  =  a^,  X.(9,o)  =  a^  ,  b(o)  =  b(o)  =  0 

(27)  \      A^cosG  sin9d9  =  0 
~o 

In  (23-27)  the  bars  have  been  omitted  and  tv;o  new 
fimctlons  g(e,t)  and  h(Q,t)  have  been  defined  by  (2[|.,25). 
The  quantity  k  is  defined  in  (21).   x:.quation  (25j)  is  a 
consequ-cnce  of  both  (7)  and  (8),  and  the  symbol  P,  in  it 
denotes  the  Legendre  polynomial  ?-,(cos9)  =  cos9. 

The  ener[:y  equation  (20)  beco:;ies,  in  dimens  ionless 
variables  and  omitting  the  bars, 

(28)-  }-    (     \^   (\.+   b[  cos9  +  -^4^  A.];4(X,9,t)  sin9  d9  + 


4- 

'o 


""   "  ^-35,  ,1  r""  .3_.^n  .n^l-6_--3 


+  |(1  -CJh){      X^sin9  d9   +'^-"^^k(|  C   X^sin9  d9)^"^='oL' 


■'..:. 'J I  L  ^y 


,  .    J. '. 


J    ■■   i: 


A 


"T' 


•■>-i.  - 
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It    should   be  .noted    that    a      and   a      are    no   lonf3;er    inde- 

o       o  '^ 

pendent,  since  they  are  related  by  the  following 
equation,  obtained  by  setting  t  =  o  in  (28). 

(20')  |. 3  l^,     a3..Si-36, 3(1-5)^^-3 

Here  <pc-  is  the  constant 

(29)  Z   =  ^   1/2 

5 

V  Method  of  Solution 

In  order  to  solve  (23-27)  for  4(3", S,t),  X(9,t), 
and  b(t),  we  make  use  of  the  fact  that  a  solution  4 
of  (23)  can  be  expressed  as  a  linear  superposition  of 

(  Tn4-1  \ 

zonal   harraonics      r~  ?    (cos9)    and   that   X   can  be    ex- 

m 

pressed  as  a  linear  superposition  of  Legendre  polynoralals 

P  (cos9).   The  coefficients  in  these  expansions  are 
m 

functions  of  t.   Next  vje  inake  the  assumption  that  all 
these  coefficients,  as  v/oll  as  b(t),  can  be  expressed 
as  series  (convergent  of  asyrii^itotlc )  of  poi/ers  of  the 
parameter  c^.   This  parameter,  defined  in(21),  is  propor- 
tional to  the  acceleration  of  gravity  g,  and  is  equal 
to  the  ratio  of  the  unit  of  length  L  to  the  hydrostatic 
head  P  /^   at  the  initial  bubble  center.   Thus  v.'hen 
0~=  0  the  present  solution  is  just  that  of  the  splrierical 
bubble  theory  since  cJ'=  0  is  equivalent  to  g  =  0. 
Formally  vie  assume 

CO  OO         t   ,     ^  /    ,  -1  \ 

(30)  4(r,9,t)  =  r  S^'Z.    ^nm^^'  r"^"'^^  ^P^(  cos  9) 


CD         OD 

..n  V- 


(31)  X(9,t)  =>   ^-"  ^>    a   (t)  P  (cos9) 
n=o     m=o 

00 

(32)   b(t)  =  y  ,^^  b  (t) 


'■-    n' 


n=o 


9. 

The  problem  is  now  reduced  to  that  of  determining 

c   ,  a   and  b   which  arc  all  functions  of  t ,  in  such 
nm'   nrn      n 

a  way  that  (30-32)  satisfy  (2[!.-27),  since  (23)  is  already 
satisfied.   Insertion  of  (31,32)  into  (26)  yields  (dot 
denoting  t  derivative) 


(33)  a   fo)=a,   a   (o)=a 
^^'>'       oo  •  '     o'    oo^       o 


a   (o)=o,    a   (o)=o    if  n+ra>o 
nm^       *    nm 


(3-';)    b^(o)  =  0,   b^(o)  =0,   n  >  0 


V;e  now  insert  (30-32)  into  (2ij.,25)  and  equate  coefficients 

of  's"^   P  (cos9),  thus  obtaining  a  sequence  of  equations 

rn. 

involving  the  coefficients  and  their  first  time  deriva- 
tives . 

The  zero  order  terms  (n  =  o)  correspond  to  th:.  solu- 
tion obtained  ".jhen  C-   o,    i.e.  in  the  absjnc^  of  gravity. 
The  equations  for  these  coefficients,  with  the  corros- 
ponding  initial  conditions  from  (33)  can  bo  solved  at 
once  and  "''■iold 


(35)  Cq.^  z   o,   a^.^  =  o,   for   m  >  1 

b   =  o 

o 

The  other  tv;o  zero  order  coefficients  satisfy 

2   • 

(36)  C   =  -  (a   a   ) 

oo       oo   oo 

(37)  (a^   a   )-a^  -  ia   -  ^  a   -  1 
^•^     ^  oo   oo '   oo    2  oo   f-  oo 

±c 

Equation  (37)  is  the  v;ell'  kno^m  equation  of  the 
classical  theory  of  spherical  bubbles,  expressed  in  terms 
of  our  units,  r.nd  (36)  is  the  corr os^oonding  expression 
for  the  potential.   The  solution  of  (37)  satisfying  the 
initial  conditions  (33)  Is  a  periodic  function  of  t 


■,y 


\ , 
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v;hich    increases   raonotonlcally   fro'n  a   rainimum  value    to 
a  maximu";!  value   during  a   half  period  and  decreases 
raonotcnically   to   the   minimum  during   the   next   half  period. 
The   solution  can  be   v/ritten   in  the   form 

(38)        a^^(t)   =  S  '\(]fy.t   +r^     =  a(tj), 

where  the  phase  t   is  defined  by  o-  ociC    )    =   d 

The  function  ^{V)    is  an  even  periodic  function 
of  '■-    defined  over  the  first  half  period  by  its  inverse 
function  X{<7C) 


(39)  ^       f  dx 

T(a)  =  V  ——-——-———    2^-.<^<<^ 


_^,s]  x-^-l-kx-^^ 

The   ntxmbers  oc    ando<  ,oc  <  oL^    are   tho   minimum   and  maximum 
values    of  «>',    and  aro    the    tvjo   roots   of    the   equation 

(i;0)      1    -   x^    -   kx"^^^-^^    =   0 

The  pararaeter   k,    xvhich  dotcr.i  nos^^and'^,    is 
defined    in    (21)    and    satisfies    the    inequality 

(UD      0   <  k   ^  k   =  -i-  (^^] 

5-1  'x  ay 

The  equality  occurs  only  when  a   =0  and  P  -   kf  •%-A'^  ,' 

The  second  cijndition  means  that  the  initial  bubble 
pressure  equals  the  iiiitial  hydrostatic  pressure  at  its 
center.   From  (l+O)  we  also  have  the  inequality  [oLis 
defined  xn  ( 29 ) J . 

(!|2)   0  <<?C(k)  <x<.-I(k)  <  1 

The  equalities  holds  \Aien   k  =  k  in  which  case  we  have 

/\ 

('4-3)   o(,iV)    =<oL=   constant 


11. 

Then  from  (38)  a^  (t),  viiich  we  will  henceforth  call 
a(t),  is  constant  and  equal  to  one.   Thus  the  unit  of 
length  L  is  the  equilibrium  bubble  radius,  i.e.  the 
radius  at  which  the  bubble  pressure  would  equal  the 
hydrostatic  pressure  corresponding  to  its  initial 
center,  while  the  total  energy  E  vjould  be  unchanged. 
In  order  to  lar-ite  the  equations  for  the  hi,'^her 
order  coefficients  (n  ^  1)  it  is  convenient  to  introduce 
the  functions  g(9,t)  and  h(9,t)  defined  by 


^  CD 

"a  mi 


ikk)     g(e,t)  =   g(e,t)  +  >_  cr^'Cb  Pi  ''^-    f  ^nV+  2i 

n=l     ^^  ^        m=o     ^''^     '     " 


+  (m  +  1)  a-^^  -"   2)c  ]   ?    )    =Z  ^Z      S  P 

nm   m  ~   °nm  xa 

n=o   m=o 


(1+5)  ii(e,t)  =  h(9,t)  -n  cr'^i   -35  i  a"^^-^  a^^  + 

n=l         k         ^° 


+  y     [i^  a   +  a  i   +  a-^-^ll:-   ]  P  )  = 
■^—_  d  nm      nm  nm   m 

m=o   a 

00  CO 

rFo      m=o  """^  ''^ 

The  functions  3(©,t),  h(0,t)  cccurrL\g  in  (iiij.,i|5)  '^re 

defined  in  (2li,25).   The  coefficients  g   ,  h   arc  defined 

nm   nm 

by  the  above  equations.   It  is  easily  verified  that 

g   and  h   depend  only  upon  these  coefficients  a   , 

^   ,  b  with  s  <  n. 
sm'   s 

Wow   inserting    (30-32)    into    (2i|.)    xve    obtain   for 

n  >   0 

e        =  -  1_  a^^'^"'^    [a        +   2  -^  a        -   g      ]      ,    m  :f  1 

nm  ,,  nrii  a     nm       ^nm        '         ' 

r,i+l 


(i+6) 

r        _ 
'nl  "2 


^   -    =   -ia3[   b^   .  ^,1   ^  2|a^^i   -   g^,]    ,   m  =  1 


Similarly    inserting    (30-32)    into    (25)    and  making  use    of 


:■'  .J 


;t      ■. 


f     =       ■■] 


•.  ;    r-  ■''  '■ 


y    J 


u...  •    ..'^\ 
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{L\.S,l4.6)    we   finally  obtain  after   so:nG  raanipulation 


•   • 


(i+7) 


aa        +  3aa        +    (a  +  35  J  a~^""-'-)a        =  ac;        + 
no  no  <;  no  ^nc 

{^\)'    =    (a3[g^^   -   a^^J)    4-   2a2h^^^      ,  ra  =  1 

aa        +  3ia        +    (1-^)    aa        =  ag       +    (m+2)ag 
nm       -^     nm  nin         °nn        ^  ^nm 

+    (ra+l)h  ,  i-n>l 

niii 

Equatiuns  (I4.7)  and  the  initial  conditions  (33,3^) 

dcterpiine  the  a   (rri4l)  and  b   once  a,   and  b,  (k<n)  are 
nra   '        n       l-.n      k 

knovm.   The  a  ,  arc  determined  in  terins  of  the  a,   with 
nl  kr,i 

k<n  by  (27).   The  C^^  arc  given  by  (1|6). 

Now  all  the  a   (niil)  arc  zero,  and  b  =0,  by  (35), 

and  a  ^  =  a(t)  is  given  by  (38,39).   Thus  wc  ms.j   proceed 

to  find  the  a^   and  b^  and  continue  step  by  step  to 
Ira      1  1^1 

obtain  all  the  coefficients.   It  v;ould  seora  from  {I4.7 ) 
that  at  each  step  n>o  uc   have  an  infinite  sot  of  non- 
horiogeneous  ordinary  differential  equations  to  solve. 
However  it  is  shown  in  Appendix  I  that  a  ,  ,  g   r.nd  h 
are  zero  whenever  m>n  or  m+n  =  odd  integer.   Thus  those 
of  equations  {I4.7)    in  i>iiich  these  quantities  appear 
become  homogeneous  linear  equations,  and  since  the  initial 
data  arc  zero  by  (33,31]-)  the  solutions  of  these  equations 
are  zero.   Therefore  we  have 

(i^.8)    a  .  =0    if  m>n  or  n+m  =  odd  integer 

b    =0    if  n  =  even  integer 
n  '■' 

Thus  from  (1|6) 

(U9)        C        =0         if  ra>n   or  n+m  =   odd    inte.'Tcr 
nm 

Hence   at   each   step    there    is    only    a   finite    set    of   equatioiis 
(exactly    [  ^^^    ])    to   be   solved. 

As   a  consequence   the    oxp??Jisions   for  4,    ^,    ^^^  ti. 


,    I 


;  !    .  1  ,  ■ .       .") 
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up    to    the    third  power    of  cT,    siiiiplify    to 

(50)  4   =  -a^ar"^   +  (^(e^^r'^P^)    ^ 6^'- 2cf'^^Z2^''^^ Z^ 

^T^(^31^'^^''  33^"^"'3^    ^   •    - 

b   =  c^^   +  (j^bo    +    .  .  . 

3 

In  the   expansion  for   \  the    terms  (^'a,,?,    andcT  ao-,P-, 

were  omitted,  since  by  (27)  a,,  =  ^-^i-  *-*•   •'-"^  ^^  "^^^ 

true,  however,  that  a_,  =  0  for  all  n.   In  fact,  for 

-"'27  ~  -1 
n  =  5  we  have  a^-,  =  -  — -   a  appa^-^  . 

It  is  convenient  to  take  advantage  of  the  preceding 

2 
results,  at  least  up  to  terms  m  cT  ,  in  order  to  calculate 

g ( 9 , t )  and  h ( 9 , t ) ,  whi ch  t hen  be c one 

(51)  i(9,t)  f  -  3a-\(a2o-  \   a^^)  '^^ {^^    '   ^  -'^^22^3  ^  -*■  • ' 


h(e,t)  =  ap^<=r+(b^  +  i  b^)^  .  I  ^2p^  2  ^ 

\    (5a2Q  +  2a22  +  ab^  [a"^a22]  *>i<^^  + 


+  \   (3a22  -  a^^[a"-'-a^2^' )  ^3^'^  +  ... 


These  expressions  facilitate  v.-riting  equations  (ii?) 
explicitly. 


-.  '^:-' 


i'  i  +  ■"■^> 
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VI      Sulutlon   Jp    toC-^ 

3r  n  =  1   the  >.'nlv  n.,.n-zci'-.  c  ..'Gf  f  iclcrit   t<    be 


» 


deter. -inod   is   b-,  .      Since   iro- :   (51)    ^■■^''^  have   ,3^-^= 
h-,,    ~   a,    Gqu:;tiori    ( !;? )    yiold^;    f.  r  b. 

(52)  (a^b^)*  =  2c? 

ThU3j  Uoin:-  the  initial  dc.tr.  in  (3'.0,  we  have 

t      o  r~,t    .        7  -^ 

(53)  "^1  =  2!:     c.-^  \  \      V?  d  X   jdT 

^o  L-'o  — ' 

equation    (53)    is    the   l-errin,;-  rite   for-nula   vhich 

acccmnts    .for   the    ■.nly   first   order   effect.      The    bubble 

shape    is    ijinaffected   in   first    er.'er,   out    the    v.r.ter 

flei-j   is   r.iodified   since 

r  r^     3 

<--j_^  =  -  ',  a^  dt. 

J 
For  n  =  2  the^only  coefficients  to  bo  detcrriiined 

arc  apQ  and  ap„.   Pro'-i  (I|.7,5l,53)  vje  have  the  equations 

/  r-i  \    "      -,  •  *    ,/■",--  k   -35 -In     _  -u   .  1  1-2  , 
(5i|)   aa^g  +  3aa2^  +  (a  +  jo  ;r.  a  ^    )a2^,   ^1   T^  "1 

a''^)=  a'^^:^  0 

(55)  aa^.  +  3^1^22  -  aa^.  =  "  ^  ^^f  >      ^   zZ   =  "^  22  "  ^ 

In  ter".i3  of  the  solutions  of  these  equations  end 
makin'--  use  ^f  (I4.6),  the  only  non-2:ero  c-.cf f icionts 
'^2in  ^"^^  given  by 

(56)  C^Q  -  -  (a^a^Q)* 

(57)  S2  =  -  ^  ^^  (a%2)* 

For  n  =  3  the  only  coefficients  to  be  deter-iinod 
are  a^^  and  b^.   From  (/;7,5l,33)  we  have 


/  .     -     \ 


iD/. 
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(58)  aa33  ^  3c.^^^  -   2ac.y^  =  -   Sb^Ca^^  +  i  it^-^a.^), 

(59)  (a^To,)*  =  -  3(a^b^a2^. )'+  6:/a2y  +  |  (a^b^a^^)*. 


b3(o)  =  b,,(rO  -  0 


Then  using  ([|.S) 


(60)  ^..  =  -  i[a^b-  -i-  3b.a^(a^,,  -  -^  a   )] 

(61)  ^^.  =  -  i[a^(a"a.J*  -  i^  b.a'^a..] 


'The  procading  equations  c"-':.  be  solved  c?cplicitlj 
in  the  sooci  1  case  of  tr.o  so-c.  Had  ''oquili'jr-.ur. 
bubble'',  v/hich  is  trc-t.-d  in  the  ne.:-,t  section.   In 
f^enerr.l,  hov;ever,  it  secns  necess  .-.r--  t  o  inteprrte 
the  preceding  equations  nunsricallv  in  order  to 
obtain  ouantitativc  results.   Such  results  ■...ill  be 
published,  soon. 

It  is  useful  to  recojnxze  that  equation  (53) 
for  the  rise  b-,(t)  cm  be  _:ait  in  the  form 

(62)   b, (t)  =  S  t^  +  P  t  +  Q 

Here    3    is   a   constant   while   }   and    '    are  periodic 

functions   of   t   having  the    3a^-.ie  pori':.>d   as   a(t).      This 

forn  ■pGx'r-.its   us   to   dctcrriine   the   behavior  of  b-,(t) 

for   all   t   by  ."icrely  cvalvu/cin;:;    integrals   for   one  period. 
These    i.atc   :i-.:l£i    appear    J.n   the    definitions   of    3,    i' 

and  "  i.hlch  are,  if  t  _  is  the  period  of  a(t) 

'"Penny  and  Fr";ce  h.'.ve  cxar.ii.ned  tlo'.  solution  of  the  hovao- 
peneous  equations  (Z^.?)  for  n  =  1  v.-ith  n(.n-scro  initial 
data  in  order  tc  study  the  p,rov;th  of  initial  asyraietry. 


16. 


(63) 


s  = 


—  t 

1  r  ^^ 

to) 
o 


u)   = 


2  ^t 


'dT 


1  { 


•o     J 


o 
/t 

\  ! 


-3 


A 


r\X-  \\    [     a-3. 


o  -/ 


^^j 


':i  =  b^(t)  -  3t^  -  Pt 


The  periodic  it"  of  ?  is  obvious,  while  that  of 

proved  with  aorie  -lanipul  tio:i.   iiqut. '-.ions  (63) 

only  xf  T  =  A   =  a  =0. 
o     o     :i 


can  oe 
aoolY 


VII   Tho  liqullibriu.!  .bubble 

In  occtiwii  V,  foli'.^KMng  o ouati •■■•:( '4.I )  it  is  shovjn 
that  if  k  =  k  then  the  zero  '..'rder  bubble  radius  a(t) 
regains  constant  and  equal  to  one.   This  occurs  if 
the  initial  radial  velocicy  is  zero  and  if  the  initial 
bubble  pressure  equals  the  pree.:.ure  in  the  v;atcr  at 
the  deptji  of  the  bubble  center,   V/c  refer  to  this 
case  as  that  of  the  equillbriurii  bubble  since  the 
bubble  vjould  remain  at  rest  in  the  absence  of  gravity. 

Since  the  radius  a(t)  is  identically  constant 
and  in  fact  equal  to  one,  due  to  our  choice  of  units, 
all  the  equations  of  the  preceding  section  siraplify  to 
such  an  extent  that  they  can  be  solved  e:rplicitly. 
By  solving  them  vje  obtain  the  following  e;:pressions 
for  X(9,t),  b(t)  and  ^{v ,9 ,t)   up  to  and  including 
the  third  order  in 


(6I4.)   X(9,t)  =  1  +(crt^)^  jf(r)  -  3A  P^vcosQ)! 


+ 


^  P^(cos9)  -I-  .  . 
P   3 


(65)  b(t)  -cr 


...2 


-  (o-t^-)    g(r)  +  .. 
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(66)        4(r,&,t)    =   -otr"2i>    (coc9)    -  ^ft   •  i^^^:^  )'^^r'^ 

12       _....  vr7::L.'i^3  _ -a  -.  ,      .^  3 


(36)^/'^ 


y3o\  icr^   r"'^  ?,  (co^G) 
;  i 


V 


In   equations    (61;,65)X=  ^i36  't   and  f(T),    gCD 
are   given  b;/ 

(67)  f CD   =  -^  -  ii^   (1  -   ccsr) 

(68)  c"(t:)    =  lo  -^  -  ^    [2(1   -    ^^^T)    -t^infi 


TheGe  functions  both  dccrov.se  fairl;/  slowly  fron 
f(o)  =  1/6,  c(f^)  =  7/1^  to  f(co)=  0,  g((£>  ,  =  9/2C  . 

Usinr   equations    (6'4.)    and    (65)    '-'e   have    cc.istruct- 
ed   graphs    of   the   bubijle  profile    at   varicua   ti.iea   froiu 
t   =  0  until   t   =  1.39,    choosing  C"  =    .222.      These 
graphs,    in  ?i~irc   I    ,    sYica-]  the   rise   of   the   bubble    as 
v/ell   as    its    deforniation. 

Since    (6li.,63')    P-re    only   the   first    tor-'is    in   a 
series    solution,    we    can   cxpocr    thern  to  roproi'.ent    the 
solution   accuxrately  for    .:    short    timo,    after   which   the 
omitted   terr.is   bocoriiO    important.       In   order   tc    determine 
v;hen   the   oiaitted  terras    becoir.e    significant,    v/e    can 
exai'iine   the    ejqoressions    in    (6!4.,6>)    to   find   when   they 
fail.      Tvjo   kinds   of   failure   cccuj?.      In   the    first   kind, 
the   oririn  of   coordinates,    v.-hlch  is    located  at   the 
center    of   ,'3ravity   of  the   bubble,    z\s.j   cross   the   bubble 
surface   and  enter  the   water.        Then  the  potential 
function  i^{r,9,t)    vjhich  has    singularities    at    the   ori- 
gin,   will   be   singular   in  the    ,:atcr,    contrary  to   the 
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assumption  of  regular  it  7  on  which  the  solution  was 
based.   ThGroforG  the  o-qTrossions  (6',l,55)  certainly 
beco':'.^  invalid  as  soon  as  X  =  0  x'ur  so  le  values  of 

0  and  t,   rr  c5>  .222  this  failure  .■.3   the  firat  to 
occur,  and  it  occur c  v;hon  (jt'~ z:  .0   ai.d  for  9  =  u. 
Thus  in  this  case  the  botto  i  of  tj;0  bubble  bulfj.os 
upv.'ard  until  xt  reaches  the  contor  of  f.ravitj. 

■fho  sec  jud  t;;-,'ije  of  fail-or'c  ocours  bocautic  (65) 
predicts  that  the  oubblc  v:ill  ri3B  for  awhile  and 
then  fall.   This  is,  of  course,  physically  unreason- 
able, and  v.'c  therefore  consider  (65)  to  fail  x^rhen 
b(t)  =  0,  i.e.  s/ren  the  bubble  stops  risin.3.   Tliis 
typo  of  failure  is  ti-c  first  to  occur  ifc:r<  .222, 
and  it  occurs  v/hen  eft  "^.8.   Therefore  v;hcn  o*  =  .222 
both  fallu_res  occur  s  i-Taltaneousl",  and  this  occurs 
at  t  =  lJ-9 ,    x-.riich  explains  -why  the  graphs  in  Figure  I 
were  only  constructed  until  this  value  of  t. 

As  an  exeiirole  of  these  results,  let  us  consider 
a  bubble  of  one  inch  radius  initially  located  ten 

inches  below  the  surface.         i:  j,- 

(A   =1   ,  Z,  =  10   ) 

o       '      o 

The  press ajT-e   above   tho    surface   is   ati-iospheric    (p      = 

1  atmosphere)    so   that   the  hydrostatic  Iiead   at   tho   bubble 
is    33    foot   plus    iC    inches    or   li.06    inches,    w>iich   is 

l\.06   tines   the   initial  radius.      Thus  cD' -  tn^o"  ~    .002:1.6. 

The  unit  of  length  L  is  equal  to  the  initial  radius 

A   for  the  eauilibriuii  bubble,  go  L  =  A   =1  inch. 
o  '  '         o 

Then  the  unit  of  tir.ie  T  =  .0025  seconds.   If  vc  use 
the  Herrinr  rise  fornula,  uhich  is  the  first  tcr;r-  of 
(65),  the  ri.:..  b(t)  =  CTt^,  so  the  bubble  i^iill  reach 
the  siu?facc  v;hcn  o  t   =  —  =  lC.;r  t  =  6ii  units. 
However  oq  aticn  (61])  falls  when  \'\r.,t)    =   0  '.•:hich 
occurs  at  t  -   16  units  or  .0'|0  secon.ds,   Tho  bubble 
rises  only  .^JS   inches  before  the  for^iulac  fail. 


'^'\:- 


'.JC'J. 


•;  :   '•  r  - 
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VIII   Disci.issii..i'.  of  rk'sults 

i:c  iaavo  'proscA''cod   a  sy^te  v.a.tic  ;.ot]ii.)d  for 
do ti^ raining  the  shape  and  rise  of  the  bi.ilblo  a;:o.  the 
potential  function  i  .^  the  vater,  frori  v/hich  the  p.r::ssure 
can  be  co..!puted.   This  Method  1g  bac-cd  "apor.  a  povjer 
series  expansion  in  torins  of  a  dlncnsionle.'.is  I'jara- 
raeter  c  which  is  the  ratio  of  the  equxlibriu"!  bubble 
radius  tf.  the  hydrostatic  head  above  tlio  initial 
bubble  center.   Thus  c"  is  prcocrtional  to  the  accel- 
eration of  gravity  g,  so  our  expansion  can  also  be 
considered  to  be  in  povjors  of  g.   For  any  bubble 
which  does  not  vont  (i.e.  break  throu-'h  the  vater- 
air  surface)  during  its  first  expansion,  c?"  if  less 
than  one  and  greater  than  zero. 

The  terns  of  zero  order  in  C   coi-respond  to  tlie 
classical  theory  of  a  sph.erical  bubble.   They  describe 
a  solution  in  ivhich  the  bubble  revialns  spherical  a:ad 
perfor".is  undajcped  periodic  radial  oscillatioixS,  •.-.Mile 
its  center  re  .ains  at  rest.   The  eqv.j.libriu  ■!  bubble 
is  a  special  case  Mliich  does  not  eve:;  oscillate. 
V..'l^en  ter/as  of  fir:>t  order  are  also  i;icludedj  the  .'nly 
inodif ieation  is  that  the  centei;  of  the  bubble  rises 
according,  to  the;  Hcrrin;;  rise  fer.iula, 

1/lien  ter:"3  of  seco.  d  order  are  incJudcd  the 
shape  of  the  bubble  is  found  to  change  althou,'j;h  the 
rise  formula  is  uncodified.   The  changes  in  the  bubble 
are  of  tv;o  kinds  corrospondi]"'g  to  tlae  two  ter'as 

2  2 

^  a2Q(t)  and  cJ  c^^it)    P^(cos!-i)  i-i    (so)  respocti^roly. 

The  first  terra  corresponds  t  o  a  change  in  volurio  of  the 

bubble,  vjhich  is  an  increase  in  the  case  of  the  equi- 

libriUi",!  bubble  and  also  in  a  numerical  example  v:hich 

v;e  have  treated.   Such  an  increase  is  to  be  c;.:pected 

in  general  on  pJiysical  grounds,  s^nce  the  bi;bble  is 

rising  into  a  region  of  lovjxr  pressure  .   Fron  the 

differential  oquation(5^)  for  a.  (t)  it  can  be  shown 

that  in  genei^al  a^^    starts  at  .^lei-o  and  increases j  at 
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least  in:..tially,  in  agrecncnt  v;itl?  the  preceding; 
cc-nsidoratic-ns.   As  a  consequence  of  the  increase  of 
a^-.(t)  it  follows  tliat  the  ti^ic  botuccn  success ivo 
rainir.ia  (or  naxina)  of  the  bub'oio  volu'.ie  is  decreased 
bolov:  that  ^^;;ivcn  h-j   the  zero  order  ( spliorlo/.'.l  bubble) 
tcr:'.:  alone.   In  addition  the  expansion  phase  is  in- 
creased In  duration,  v;hi].e  tnc  contraction  phase  is 
dirainlsh.cd  in  duration. 

The  other  second  order  tcr.'i  corrosnonds  to  a 

flattening  cf  the  bubble.   Since  a^^(t)  starts  at  zero 

22 

and   dccroasos,    o.s    can  bo   seen  iron  the   differential 
equation(  55) ,    this    flattening   is    in^  t}io   vertical 
direction,    i.e.    the    vertical    separation   botv/een   top 
and   bottor.i   is    di-r.i  -ished,    iSiile    the   horizontal 
separation  betivecn   sides    is    increased.      This    can  be 
seen   in  the   figures,    and  is    in  arreci;i.ent  with  ob- 
sorvaticn,    as    are   the  procedine  results, 

Vfncn  terivis    of  third   order    In  c  are    included 
the   shape   of   the    bubble    is    arain  .-lodified  while    the 
rise  fornula   is   also  chonf^ed.      The    only  t&r"i  of   this 
order   in   the   equation  for   the   bubble  radius    is 
(f^   ao-.(t)    ?-(coafe),    vrhich   leaves    the    bubble   volurae 
unchanp:od  but    chan^'os    its    shape.      Since   a-,^(t) 
starts    at    zuro    and    increases,    as   one    can   show  fron 
the    differential    equation    (5^-),    this    chan,;.' e   of   shape 
corresponds   to  pullf.nc  up   the   botton  and  pushing  up 
the   top,    so   the    bubble   becoines   kidney-shaped  x-;ith  the 
lovrer   side   concave.      This    behavior    can   be    seen   in   the 
flf^iires    and    is    also    in  agrceraent   v/ith  observation.       The 
third   order   correction  to    the   rise    for.ViUla   can   be   shoi-m 
to   start    at   zero  and  to  becor.ie  nerativc,    fro'':i  equation 
(59),    in   a"rcer.ient    udth  the   numerical    example   and   the 
oquilibriixu  bubble.      This    correction   is    also    in    a;f^rcc- 
uient   with  experiment,    since   the   liorrinj  forr.iula  pre- 
dicts   too   large    a  rise. 

All   of   the    gravity   effects    dosci-ibed   here    are 
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greater  in  ::iaf\:nituJc  ^ar-c  occox'     \ctv    ouicl-.-ly,  the  larger 
tho  value  of  cT.  '.'ou  ^   is  lat^'xr  fer  largo  cicplos ions, 
so  all  the  ciTccts  cU.3cribod  should  bo  'ioro  pro  -li/icnt 
for  such  oxplosio;^s.   A  l.:.r2-  vali;e  of  (.^  is  also 
obtalnec:  if  the  pressure  above  tho  \;ater  surface  is 
reduced.   This  explainn  vzhy  these  effects  have  been 
observed  m  reduced  preseure  tanks. 

The  behavior  of  tho  bubble  depends  not  only 
upon  c'  'out  also  upon  another  di^iensionless  paranctcr 
k.   Tho  s-ialler  the  value  of  ]i  the  ,';';rcator  is  the 
araplitudo  of  oscillation  of  the  bubble,  since  si-;all 
k  corresponds  to  a  larp-e  e::ploslon  energy 


K  r 


6-1  - -5  ^ 


0-1 

A 

The  raaximuiu  value  of  k  is  k  which  is  attained  for  the 
cquilibrixii'a  bubble.   Since  the  ^:;:ravitational  effects 
depend  raainly  on  cT,  rosultr,  for  tho  oqiMlibriu..; 
bubble  ir.dicate  vjhat  vrill  ho.pper.  in  reneral.   There- 
fore v/e  c:;poct  that  f..vr  any  value  of  k,  ^quatit.ns 
(50)  for  \   and  b  x-;ill  fail  at  a  certain  tine,  just  as 
they  do  for  the  equilibriuu  bubelo,  k  =  k.   Since  in 
general  the  bubble  becovaes  flattened  and  then  k^.dne:/- 
shaped,  \ju    c;rpect  the  center  part  of  i;  he  loT-Jcr  surface 
to  continue  rising,  like  a  jet,  until  it  reaches  the 
upper  sui'^face,  thus  converting  the  bubble  into  a 
torus.   For  the  equilibriu.-.!  bubble,  failure  of  the 
equations  occurs  whcncrt"~.S,  so  \:c   r-iay  assuno  that 

the  real  breakdevjn  for  any  bubble  vail  occur  when 

2  -l/? 

iTt  i:rl  or  t^j-  '".   As  one  period  of  tho  bubble  is 

about  3/2  units  for  an  er-gp lesion  bubble,  the  breakup 
should  occur  after  about  ^cr-^/2  .^^^^^^^^      p^^  ^.^^^ 

OT'losions,  with  large  values  ofcr,  this  nay  occur 
diiring  the  first  period. 

If  tho  bubble  docs  breakdovni  into  a  torus,  '.hat 
v;ill  be  its  subseqvxent  behavior?   The  results  of 
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Toractika  on   a  torus  of  ■.mc   fluid  -/itbin  another  fluid 
indicate  that  a  torus  is  also  unstable  and  v:ill 
generally  breakup  into  a  certain  n-iJt;,ibGr  of  pieces 
dctcr:".incd  by  the  dimensions,  etc.   Then  eac:\  of  these 
pieces  rairht  flatten,  becorie  kidney  shaped  aiid  break- 
doi>m  intj  another  torus,  whicJi  ijould  again  breakup, 
etc.   This  socvris  to-  1-appen  when  a  drop  of  one  fluid 
(e.g.  inic )  is  dropped  into  anv„thor  (e.f!;,  uator). 

Finally  let  us  consider  xiou  corprcssibilj  ty  will 
alter  the  pr^-sont  results.  via   l\c.vo   c.ttcr.iptcd  to 
account  for  c^r.pressibility  of  the  uatcr,  as  ivc  did 
in  our  previous  report  on  the  spiicrlcal  bubb].e,  by 
using  the  wave  equation  for  the  potential  in  the 
water  rather  ti'.an.  Laplace's  oquatiori."  '..'o  then 
expressed  the  solutioji  of  this  equation  as  an  infinite 
scries  of  rralti'^Liles,  but  v:c  wore  forced  to  use  a 
fixed  ori,"ir.  rather  than  a  noving  one.   V.'o  also 
expressed  the  bubble  surface  as  r,    scries  in  pov.;er3  of  c" 
and  in  Lcgondre  polyno:nials,  and  atteviiptod  to  dctor;:iinc 
the  coefficients.   Uo  found  that  the  terv,:.,s  of  zero 
order  in  c?"  wore  exactly  the  sa.ve  as  the  stlution  of 
our  provloiis  report,  as  we  naturally  e:qD acted.   The 
higher  order  coefficients  satisfied  second  order  equa- 
tions sir.iilar  to  those  obtained  in  the  present  report, 
but  the  coefficients  in  tlicsc  equations  involved  the 
function  a(t'  of  our  previ^ais  report,  rather  than  the 
incompresslblG  solution  a(t).   'fhis  was  not  the  only 
difference  but  sec::ied  to  be  t/ie  i.iost  i.iportaiit  one. 
vJo  therefore  G::pect  that  all  of  the  present  results 
will  still  apply  in  a  corprossibl..  fluid,  and  that 
the  v.iain  effects  of  coirrpressibility  v;oulu  be  accounted 
for  by  replacing  a(t)  in  all  equations  by  the  conpressiblc 
solution  of  our  previous  report.   In  this  v;ay  the 
naiii  result  of  that  repcrt-da?rrped  radial  oscillation- 
would  be  incorporated  into  the  present  work, 

"In  an  unpublished  vianuscript . 
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Appendix  I 

Theorem:   If  i|),  X  and  b  are  represented  by  the  formal 
series  In  equations  (30-32)  and  If  they  satisfy  the 
equations  (24-27)  then  for  all  n  >  0 

(Al)   g„„  =  h„^  =  a^^  =  c„^  =  0 

°nm    nra    nra    nm 

if  m  >  n  or  If  ra  +  n  Is  odd,  and 
(A2)   b^  =  0 

If  n  Is  zero  or  even. 

Furthermore  g   and  h   are  polynomials  In  a,  ^  , 
nm      nm  k  jj 

&,  j^  ,    b,  and  b,  with  0  <  £     <  k  <  n  and  are  Independent 

of  all  other  a's  and  b's  except  a^^. 

oo 

Proof:   The  proof  relies  on  induction  with  respect  to 

n.   Thus  we  assume  that  the  theorem  Is  true  for  all 

n  <  N  and  we  will  show  that  this  implies  the  theorem 

for  n  =  N  +  1.   It  has  already  been  shown  in  the  text 

that  the  theorem  is  true  for  n  =  0  (see  equations  (35) 

and  (51))  and  therefore  the  theorem  will  follow  for  all  n. 

In  order  to  perform  the  inductive  step,  we  Intro- 
y, 
duce  the  set  H,  '  of  all  formal  power  series  in  a   and  x 

with  coefficients  depending  upon  a  parameter  t,  which 

satisfy  the  conditions  u   (t )  =  0  if  n  <  k,  or  m  >  n  and 

nm 


n<    >^,    orm  +  n  =  odd  and  n  <    Z  .      The  quantities 

nm 
given  by 


u   (t)  are  the  coefficients  in  such  a  formal  series  u 
nm 


00        00 

(A3)   u  =  ]~  CT^  i:  X™  u„^(t). 

-i'    £■ 

We  also  define  the  subset  Yi/    C.    H.   of  power  series  in 


which  the  non-zero  coefficients  u  ,    n  <   £  ,    are  polyno- 

nm    -~ 

mlals  over  some  set  [c] . 

If  X  =  cos  e,  then  since  P  (x)  is  even  or  odd  ac- 
cording as  n  is  even  or  odd  it  follows  from  the  induction 


:»•. 


y*r.^»    r-    Si  !ji. 


V^^'.i      />      y- i':.'      [-:'■ 


0    < 


J.  .1  >' 


i.     ■- 


tf\:: 


'5)i51      L  .,;>ri 


!^     f     f 


•-:■:/■■, 


i")0 


/^     .  a  •  1  >  3 


!P)  s: 


!./.     ..■^;:^«, 


-Oi'lv -!"'^'^    .. 


■  ••■♦   .->,■ 


•■     :!■>;■: 


.    V   'f  )     .     .^' 


-  *J 


-l-vv     J  x     j;,/,o 


24. 


hypothesis  that 

(A4)   g  e.  H^,    h  £  H^. 

If  we  let  [c]  denote  the  set  of  a   ,  c   ,  b  and  their 

nm   nin   n 

time  derivatives  with  0  <  n  <  N  then  it  also  follows 
from  the  induction  hypothesis  that 

(A5)   Mr.t)  £  H^,    A  &  H^,    bx  £  hJ. 

We  first  wish  to  show  that 

(A6)   g  £  H^-'^,    h  £  H^-"^. 


In  order  to  do  this  we  first  note  the  following  proper- 
ties of  the  sets  H,  : 

(a)  Hj^    '-    Hj^.  i^        i    >  i^  '        and        k  >  k' 

(b)  u  £  ii/  ->  u  ^a^x"^  6  H^        n  <   Z 

o  nm  n  — 

(c)  u  £  E^,        V  £   Hj^^   ->   (ecu  +  pv)   £  Vl^ ,        uv   £  h£"^^ 

(d)  u  S  Hjf ,        V  S  H^";',         /  +  k'    <   i'    +  k  -»  uv  £  iij^Jf 


(e)  u  £  H^,        Oc    ?^  0  ->   (a.    +  u)^  -  nOc""^u  £  H^'^'*'-^, 
(CX  +  u)"  -  CX"   £  H^^ 

(f )  u  £  h/-  ->  u^  £   h/' 

(g)  u  6  H^^,        V  £  H^^  -^    (1    -  x2)u^v^  £   h/""^  . 

Now  making  use  of  these  properties  and  the  induction  hypo- 
thesis we  can  prove  (A6). 
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To  prove  (A6)  we  consider  various  parts  of  g  and  h 
separately.   Thus  If  we  consider  the  following  terras  In 
g  we  have 

(A7)   ■|.^(  X,e,t)  *   f  a"  Z  [2  f  a^  *  (m  *  1  )a '"^^'o^lP^ 

n=l    m=0 

=  a^a[  A-2  +  2a~^(  A  -  a)] 


00        00 


"  ^  (m  ^  l)c   (a-("^-^2)  _  ;^-(m+2))p  ^  gN-^l 


-,  n  nm  ra    . 

n=l    m=0 

This  is  proved  by  noting  that  since  a  7^  0,  the  first  term 

above  is  In  H^    according  to  (e).   Furthermore  each  term 

in  the  sum  is  In  H*:^^  by  (e),  (d)  and  (a),  and  therefore 

in  H^"^-^.   Now  by  (c)  the  s\am  Is  also  In  H^"^-"-.   Thus  the 

statement  (A7)  is  proved.   In  a  similar  way  the  remaining 

terms  of  g,  and  all  the  terms  in  h  may  be  shown  to  belong 

to  H^"*"-^,  but  the  details  will  be  omitted.   Therefore  (A6) 

may  be  considered  proved. 

Now  we  consider  equations  (47)  for  the  a^^  and  b„ 

^  nm      n 

with  n  <  N  +  1.   Whenever  ra>norn+ra=  odd,  the  Inhomo- 

geneous  terms  in  the  equation  for  a   vanish,  since  these 

terms  are  proportional  to  h   and  g   ,  which  have  lust 
^   ^  nm     °nm  '' 

been  shown  to  vanish.  The  resulting  equations  are  linear 

and  homogeneous,  and  by  (33)  the  initial  data  are  zero. 

Therefore  it  follows  that  the  solutions  a   =0  for  m  >  n 

nm 

or  n  +  m  odd.   Similarly  from  equation  (47)  for  b  we 

find  that  b   +  a  ,  =  0  for  n  even.   In  order  to  conclude 
n    nl 

that  b   =0  for  n  even  we  must  prove  that  a  t  =  0  for  n 
n  nl 

even.   Once  this  is  proved,  it  will  also  follow  from  (46) 

that  c    =0  for  m  >  n  or  n  +  m  odd  and  the  theorem  will 
nm 

have  been  proved. 

Thus  we  must  show  that  for  N  odd,  a^  ,  ,  =  0.   To 
this  end  we  consider  equation  (27)  and  expand  the  inte- 
grand in  powers  of  a,  using  (31)  for  A.   The  coefficient 
of  a    in  this  expansion  is  a  sum  of  terms  of  the  form 
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pi  4 

^.  i=l 


(A8)  ^    f^n  m  ^m  ^^^^^  ^^• 

^,     i=l     ^i"^i  ""i 

■■^  21  n^,  =N+1 

In  this  equation  the  variable  x  =  cos  6  has  been  intro- 
duced.  If  all  n.  in  a  particular  term  with  a  non-vanishing 
integrand  are  <  N,  then  since  their  sum  =  N  +  1  =  even, 
the  sura  of  the  corresponding  m,  is  also  even.   This  fol- 
lows from  the  induction  hypothesis,  since  for  n  <  N, 

4 
n.  +  ra.  =  even  if  a_  _  ^  0.      Consequently  IT    P^  (x)  is 
i    1  ^i  1  i=l    i 

an  even  polynomial  in  x  and  therefore  the  integral  vanishes, 

Hence  the  only  terms  not  vanishing  identically  are 

those  in  which  one  n,  say  n.  =  N  +  1  and  the  other  n^  are 

zero.   Now  if  n.  =0  then  m.  =  0  since  otherwise  a^^   =  0. 
^   i  i  om^ 

Therefore  IT    P   (x)  =  P   (x)  and  the  integral  vanishes 

i=l  "^i       "^j 
xonless  m.  =  1  due  to  the  orthogonality  of  Legendre  poly- 
nomials. Consequently  the  only  non-vanishing  terms  are 
of  the  form 

(A9)         a^J   ajj^3^^3^P^(x)x  dx  =  |  a%^3^^^. 

-1 
There  are  four  such  terms  and   thus  the  coefficient  of 
(T    is  ^  *^N+1  1*   ^^°"^  equation  (27)  this  coefficient 
must  vanish  and  therefore 

(AlO)  ^^^^^^   =  0. 

This  completes  the  proof  of  the  theorem. 
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Figure  1 
Rise  and  Change  of  Shape  of  "Equilibrium"  Bubble 


Profiles  are  shown  for  five  values  of  t  from  t  =  0  until 

t  =  1.89,  at  which  time  breakdown  occurs.   The  center  of 

gravity  is  also  shown  for  each  value  of  t  by  a  dash  on 
the  z-axis. 
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